Role of line-of-sight cosmic ray interactions in forming the spectra 
of distant blazars in TeV gamma rays and high-energy neutrinos 
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ABSTRACT 

Active galactic nuclei (AGN) can produce both gamma rays and cosmic rays. 
The observed high-energy gamma-ray signals from distant blazars may be dom- 
inated by secondary gamma rays produced along the line of sight by the in- 
teractions of cosmic-ray protons with background photons. This explains the 
surprisingly low attenuation observed for distant blazars, because the produc- 
tion of secondary gamma rays occurs, on average, much closer to Earth than the 
distance to the source. Thus the observed spectrum in the TeV range does not 
depend on the intrinsic gamma-ray spectrum, while it depends on the output of 
the source in cosmic rays. We apply this hypothesis to a number of sources and, 
in every case, we obtain an excellent fit, strengthening the interpretation of the 
observed spectra as being due to secondary gamma rays. We explore the ram- 
ifications of this interpretation for limits on the extragalactic background light 
and for the production of cosmic rays in AGN. We also make predictions for the 
neutrino signals, which can help probe the acceleration of cosmic rays in AGN. 

Subject headings: gamma rays, cosmic rays, active galaxies 



1. Introduction 

Active galactic nuclei (AGN) are believed to produce both gamma rays and cosmic 
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Acciari et al.l 120091 ; iCostamante et al.l 120081 ) . It is believed that blazars also accelerate 
high-energy cosmic rays, but it is more difficult, if at all possible, to associate the ar- 
rival directions of cos mic rays with extragalactic sources because o f the effects of the Milky 



Way magnetic fields (IHarari et al. 



2000a 



b|; lAbraham et al.l 120071 . 120081 : iGolup et al.l 12009 



Murase fc Beacom 201ol ; Calvez et al. 201ol ). which are known to have magnitudes as large 



as several microgauss. Gamma -ray spectra o 



universal photon backgrounds (ptecker et al 



: distant sources haye been used as a probe of 



1992 



de Jager et al. 



1994 



Salamon fc Steckei 



19981 : Istanev fc Franceschinill 19981 lAharonian et al.ll2006HMazin fc Raue!l2007 ; iFinke fc Razzaqud 
20091 



Yang fc Wan3l2010l : 



Abdo et al. 



2010h . However, it was recently pointed out that the 



highest-energy gamma rays from the more distant objects may be secondary gamma rays pro- 
duced relatively close to Earth in proton-photon interactions, and n ot primary gamma rays 
emitted at the source fessev fc Kusenkol l20ld : fessev et alihoiObh . Although th e energy 
attenuation length of protons below the Greisen-Zatsepin-Kuzmin (GZK) cutoff (IGreisen 
19661 ; IZatsepin fc Kuzminlll966l ) is much larger than the horizon size, the protons do rarely 
interact with the background photons and produce secondary gamma r ays at a rate that 
resul t s in observable flux and shows an excellent agreement with the data (lEssey fc Kusenko 
20101 ; lEssey et al.ll2010bl ). The contribution of cosmic rays to the observed spectra of distant 
blazars should not be neglected in analyses of gamma-ray opacity of the universe and in 
deriving constraints on the photon backgrounds. 

The identification of observed gamma rays with secondary showers along the line of 
sight reconciles the observed TeV spectra with theoretical models. The primary gamma 
rays should exhibit clear signatures of absorption due to their interactions with extragalac- 
tic background light (EBL), suc h as a sharp cutof f at energies around 1 TeV that was 



predicted prior to obse rvations ( Stecker et al. 1992h. However, the observed spectra d o 



not show such features (lAharonian et al.l^Ooll lAcciari et al. 



2009 



ICostamante et al. 



2008). 



This can be explained by the lower levels of EBL (lAharonian et al.l l2006t iMazin fc Raue 
2007 ; Finke fc Razzaauel 2009 ). combined with much harder in trinsic spectra of distant 



zuim umKe cz rtazzaq uej izulm ). comDmea witn mucn naraer in trinsic sp 
blaza rs dstecker et al.ll2007h th an those predicted by earlier models (|Peacockll 



19871 : iHeavens fc Drurvlll988l : iBednarz fc Ostrowskilll998l : Malkov fc Drurv 



981 



Kirk fc Schneider 



200l|). Alterna- 



tively, the lack of absorption features can be explained by the production of seco ndary gamma 
rays, which replace the pr imary gamma-rays lost to their interactions with EBL (lEssey &: Kusenko 
20101 ; lEssey et al.ll2010bl ). We will present additional evidence that the latter approach pro- 
vides a consistent explanation of all the present data. 

Secondary gamma rays produce point images as long as the magnetic deflections of 
protons are sufficiently small. The magnetic fields in the source host galaxy introduce only a 
negligible image broadening, at most, of the order of the galaxy size divided by the distance 
to the source. In intergalactic space, the magnetic fields are much weaker than inside a 
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galaxy, and can ha ve the field strenghths as low as a femtogauss (lAndo fc Kusenkol 12010 



Kandus et al.ll2010l ). Therefore, the proton deflections are small, and the images of sources in 
secondary photons appear pointlike. Indeed, the deflections of protons emitted by a distant 
blazar on their passage in a random magnetic field below 10 fG are smaller than the angular 
resolution of atmospheric Cherenkov telescopes (ACT): 

a$ ~ o.i° i * ^ r 4 x io : g ^ ("a 1/2 r ju 1/2 . ri) 



io- 14 g; v e j ViGpcy V 1m p c 

Here D is the distance to the source D and l c is the average correlation length. 

It is easy to understand qualitatively why, for a sufficiently distant source, secondary 
gamma-ray flux should dominate over primary gamma-ray flux. The flux of the primary 
gamma rays is attenuated by their interactions with EBL and, for a source at distance d, 
the flux of unattenuated high-energy gamma rays is 

Fp limaiy!l (d) oc — exp{-d/A 7 }, (2) 
where A 7 is the attenuation length due to the interactio ns with EBL. Cosmic rays emitted 



from the same source with energies below the GZK cutoff (lGreisenlll966l ; IZatsepin fc Kuzmin 



19661 ) of about 3 x 10 10 GeV can cross cosmological distances without a significant energy 
loss. Their flux at distance d is Fprotons oc 1/d 2 . Although the photon background is optically 
thin for cosmic rays, the protons do rarely interact with the cosmic background photons and 
produce gamma rays. 

Let us consider an isotropic source of protons. Since the proton flux is not attenuated, 
the same number of protons pass through every spherical shell at any radius r. Let protons 
produce gamma rays at the rate p per unit length of path, and the number of gamma 
rays passing through a shell at distance r is $ 7 (r). Some gamma rays are lost to pair 
production characterized by attenuation length A 7 . The change in the number of photons 
$ 7 that occurs between r and r + dr is d$ 7 = pdr — (1/A 7 )$ 7 dr. The solution of this 
equation with a boundary condition $ 7 (0) = 0, appropriate for secondary gamma rays, is 
$ 7 (r) = pA 7 [1 — exp(— r/A 7 )]. For the flux per unit area F S econdary,7(^) = $ 7 (c?)/(47T(i 2 ), this 
solution gives 



-^secondary ,7(^0 Altd? ^ ^ ^ 7 ] ^ ^ 



1/d, for d <C A 7 , , . 



1/d 2 , forrf>A 



7' 



This derivation extends to beamed sources, as long as the effects of the beam broaden- 
ing a re small, which is t he case for IGMF below 10 fG, in agreement with observational 



data (IKandus et al.ll2010l ). 
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It is clear from Eqns. (J2H3J) that, for a sufficiently high proton flux, secondary gamma 
rays should dominate the spectra of very distant sources above E ~TeV, because their flux 
suppression is less severe at large distances. 

Secondary neutrinos also show a different scaling wi th distance. This c an be used as a 
tool to distinguish between primary neutrino s from AGN (IStecker et al.lll991l) and secondary 
neutrinos produced along the line of sight (jEssey et al.ll2010bl ). For neutrinos, there is no 
absorption, and the flux of secondary neutrinos scales as 



-^secondary,;; (fl) (-fprotons X oT) OC . 



(4) 



The 1/d scaling applies as long as the intergalactic magnetic fields (IGMF) are sufficiently 
small to allow the protons to remain within the angular resolution of the detector. This 
unique scaling law can be exploited by future experiments: a larger number of sources 
should be within the field of view than one would predict based on the primary gamma-ray 
flux. 

The secondary gamma rays should arrive at random times and show no temporal vari- 
ability, and this is consistent with the data. While variability has been observed for nearby 



TeV blazars at TeV energies (IKrennrich et al.ll2002t iBlazejowski et al. 



20061 ) and for distant TeV blazars at energies above 200 GeV (jAcciari et al. 



2005 



Aharonian et al. 



2009 



20101), no 



variability has been observed for distant TeV blazars at Te V or higher energies. The flux of 
gamma rays with E > 200 GeV is dominated by the photons with energies E w 200 GeV, 
not TeV. We expect to see no variability above TeV for blazars with z > 0.1 considered here. 

The current data are consistent with the interpretation that secon dary gamma rays 
dominate observed signa ls from distant sources at the highest energies (jEssey fc Kusenko 
2010l : lEssev et al.lhoiObh . 



Of course, the relations in Eqns. ([2]|3|) are crude approximations, and one must take 
into account a complex chain of interactions that ensues when the first interaction triggers 
an electromagnetic cascade. To this end, we will present the results of numerical Monte 
Carlo simulations that account for all the relevant interactions and we will compare theo- 
retical predictions to the data. The overall normalization is uncertain, but the shape of the 
spectrum is fixed. Hence, our model predictions should be viewed as one-parameter fits to 
the data. In all the cases we find an excellent agreement between the shape of the spec- 
trum predicted by the model and the data. We will discuss the implications of our results for 
studies of extragalactic background light and other areas where it is important to distinguish 
between primary and secondary gamma rays. In what follows we will assume that IGMFs 
are s maller than 10 fG; th is assumption is in agreement with all the present observational 
data ( IKandus et al.ll2010f ). and it allows us to explore the role of line-of-sight cosmic-ray 
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interactions in the case where they dominate, while the highest-energy intrinsic gamma rays 
ori ginating at the sour ce are filtered out. The effects of larger magnetic fields were discussed 
bv lEssev et aP feoiOah . 



2. Gamma rays from distant blazars 



To obtain the spectra that can be compared with the data, we employed a detailed 
Monte-Carlo simulation, tracking individual cosmic ray protons and all secondary particles 
from the source through the intergalactic medium. The intergalactic medium was modeled 
by a background photon field, consisting of the CMB and EBL, and IGMF. 

We have considered a wide range of EBL models, including the model with t he lowest 
level of background light based on t he lower limits obtai ned from galaxy counts (IXu et al. 
20011 ) and the highest level based on IStecker et al.l ( 120061 ). The IGMF was modeled by 3D 
cubes of the size of a characteristic correlation length, with the direction randomly chosen 
for each cube. 

The dominant interactions for the protons are proton pair production (PPP) and pion 
photoproduction: 



V + lb -> V + e + e " 

V + lb -> U + 7T + 
P + lb ->■ V + 7T° 



(5) 
(6) 
(7) 



where 7& is a background photon. Neutrons and pions subsequently decay and produce 
neutrinos, photons, electrons and positrons. Individual proton interactions and EM cas- 
cades were modeled using a standard Monte Carlo approach where energies and direc- 
tions were sampled from distribution functions constructed from t he appropriate cross sec- 
tions ( Szabo fc Protheroe 1994 ; Blumenthal 1970 ; Protheroe 19861 ). The outgoin g distribu- 



tion f unctions for pion photoproduction were generated using the SOPHIA package (IMiicke et al. 

2oooh . 



For each iteration, particles were propagated a distance far less than the average corre- 
lation length of the magnetic field, to ensure the accuracy in calculated deflections. Two cuts 
were applied to the particles arriving at the z = surface to decide whether or not to include 
them in the observed spectrum. First, the particle must point back to a point in the sky that 
is within an angular distance defined by the point spr ead function (PSF) of the observing in- 
strum ent. For energies below 100 GeV the Fermi PSF (IRando fc for the Fermi LAT Collaboration 
20091 ) was used, and for energies above 100 GeV a PSF for a typical atmospheric Cherenkov 
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telescope such as HESS, MAGIC, or VERITAS ([Holder et al.ll2008l ) was used. Second, the 
particle must arrive within a cone that is characterized by the jet opening angle for the 
source. 

The results for the spectra are presented in Fig. 1. We have chosen three most distant 
blazars observed in the TeV energy range, which show no variability, and which have upper 
limits set by Fermi: 1ES 0229+200, 1ES 0347-121, and 1ES 1101-232. We fit the spectra 
with secondary gamma rays produced by cosmic-ray interactions along the line of sight. 

There are several uncertain parameters relevant to this fit. For the IGMF we chose the 
average value of B = 1 fG, coherent over a correlation length ~ 1 Mpc; this is consistent with 
the observations (lAndo &: Kusenkoll2010l ). (Effects of the magnetic fields on the goodness of 



fit will be explored elsewhere.) The jet opening for the protons was assumed to be 9j et = 6°, 
which corresponds to a moderate Lorentz factor of T = 10. The choice of 9- ]<st will only 
have an effect on the spectrum if the deflections of the secondaries are greater than #j et . 
For the magnetic fields considered in this paper this is only true for gamma rays in the 
GeV range. Thus the fits to the TeV energy range are insensitive to the choice of 9- ]et . 
However, the choice of 9 )et does affect the overall luminosity normalization, as shown in 
Fig. 5 and d iscussed below. The spectrum of protons emitted by the AGN is unknown, but 
as shown by lEssey et al.l (l2010bl ). the spectra of secondary photons are not sensitive to the 
spectral index a and the maximal ener gy E max of the protons in a broad range of parameters 
consistent with the data on UHECR (IBerezinsky et al.ll2006l ). We parameterize the proton 
spectrum by a simple power law: 



F(E)cxE- a 9(E max -E). 

We allow a to vary between 2 and 2.7 because a = 2 and a 
cosm ic ray data at low energy and at the highest energies, respectively (IBerezinsky et al. 



2.7+£;?5 give a good fit to the 



20061 ). In fact, a satisfactory fit to the data can be obtained for any value of 1.5 < a < 3. 
This should be contrasted with fitting the data in the absence of cosmic ray contribution, in 
which case very hard in trinsic gamma-ray s pectra of distant blazars are required (cf. Fig. [2] 
and the recent paper of lEssey et al.l (j2010al )). 



What is also uncertain is the power produced by AGN in cosmic rays, L p . We fit the 
photon data by choosing this unknown parameter. The predictions shown in Fig. 1 are not 
sensitive to most of the parameters listed above, except for L p . Furthermore, while the 
spectra depend on the model of EBL for a fixed value of L p , the uncertainty in EBL is much 
smaller than the uncertainty in L p , so, to first approximation, it could be absorbed by the 
uncertainty in L p . For blazars, relativistic beaming effects must also be taken into account 
when calculating L p . The isotropic (unbeamed) equivalent luminosity L p ^ iso can be much 
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Fig. 1. — Comparison of the predicted spectra with the HESS data for th ree blazars: pan- 



els (a) and (b) show model prediction and the data for 1ES 0229+200 (jAharonian et al. 
2007bT); panels (c) and (d)) show the predicted spectrum and the data for 1ES 0347- 



121 ( jAharonian et al.l l2007al): pane 



for 1ES 1101-232 (lAharonian et al 



(e) and (f) show the model prediction and the data 
2007 cJ). The Fermi upper limits shown at lower energy 
were derived from the data by iNeronov fe Vovkl (120101 ) . Panels on the left show the pre- 
diction for "high" EBL, while panel s on the right show the prediction for the "low" EBL. 
The "high" EBL is from the model of IStecker et al.l ( 120061 ) . while the "low" EBL is the result 
of scaling down of "high" EBL to the level of 40%. (This range encompasses all published 
models.) 
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Fig. 2.— Same as in Fig. 1(b), with the addition of a gamma-ray signal expected in 
the absence of line-of-sight cosmic-ray contribution, shown by the short-dashed line. The 
intrinsic gamma-ray spectrum was chosen to be a single power law with spectral index 2. For 
this illustrative example, the predicted signal at lower energies violates the Fermi bound, but 
a harder spectrum and/or deviations from a single power law can improve the fit to the data. 
While an acceptable fit can be obtained without a cosmic-ray contribution, the secondary 
gamma rays from cosmic rays fit the data for a broad range of parameters without tuning. 
A broad rang e of spectral indices and magnetic fields, with and without cosmic rays, was 
considered by lEssey et al.l (j2010al ) . 
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greater than L p : 



L 



p,iso 



1 — cos 9\ 



JCt 



L 



V 



(9) 



As discussed below, a more detailed goodness-of-fit analysis does show some preference for 
some models of EBL. To encompass all t he published models using a simple parameterization, 



we have taken the "high" EBL model of IStecker et al.l (120061 ) as the upper limit and an EBL 



of models fall in this range ( 


Sa 


amon & Steckei 


1998: 


Kneiske et a 


.112002. 


2004; 


Stcckcr et al. 


2007; 


Franceschini et al. 


2008; 


Horiuchi et al. z 


>009; 


Primack et al. 


2009 


■ 


Gilmore et al. 


2009; 


Razzaque et al. 


2009; 


Finke et al. 


2010). 



The model predictions presented in Fig. 1 are essentially one-parameter fits, with the 
shape of the spectrum fixed by the shape of EBL, and the overall height proportional to the 
product of L p and the level of EBL. The parameters used for the spectra shown in Fig. 1 
and the goodness of each fit are shown in Table 1. 

We have obtained an excellent fit for each of the three blazars using so me reasonable val 



uesof cosmic ray power, consistent with theoretical models of cosmic rays (IBerezinsky et al. 



20061 ). Furthermore, we have obtained a reasonably good fit for any model of EBL, al- 
though the power in cosmic rays required in each case depends on the level of EBL. If 
anything, "high" EBL models are favored, but it will take more data to achieve a statisti- 
cally significant discrimination between different EBL models. This conclusion is different 



;he conclusions of other authors, who considered only the primary photons to fit the 



from 

data (lAharonian et al. 
20091 



Abdo et al 



20061 ; Mazin fc Raudl2007l ; lAharonian et al.ll2007bl : iFinke fc Razzaque 



20101 ) . It is easy to understand the origin of this discrepancy. If one tries to 
fit the data with primary photons alone, distant sources should show much more significant 
degree of absorption at high energies. The apparent lack of absorption can be compensated 
by lower EBL and harder injection spectra. However, for secondary photons, the absorption 
is much less dramatic because a large number of observed photons are produced relatively 
close to Earth. Furthermore, higher EBL density actually increases the production of gamma 
rays in pj interactions. 

The effect of different E max and a is to change the relative contribution of reactions 
flSJ) and (j7J to the flux of secondary gamma rays. If the proton spectrum extends to very 
high energies (as suggested by data on ultrahigh-energy cosmic rays), then the PPP reaction 
on CMB ((Sj) dominates over pion photoproduction on EBL (j71). If, however, E max is small 
(but larger than the threshold of 10 7 GeV), the pion photoproduction on EBL (J7j) is the 
dominant source of gamma rays. The only difference for the gamma-ray spectra is the power 
L p required from a given source. Since this power is unknown, one can obtain equally good 
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fits to the gamma-ray data for different values of E max , as shown in Figs. 3 and 4. The only 
differences in the photon spectra occur at energies well above the reach of data. 

Of course, the luminosity of the source in protons required to fit the observed flux 
depends on both E max and a. This dependence is shown in Fig. 5. 



3. Neutrinos from distant blazars 

Neutrino spectra are much more sensitive to the values of E mSuX and a. This is because 
the shape of the gamma ray spectrum at the end of the electromagnetic cascade is determined 
mainly by the spectrum shape of the EBL, regardless of whether the main contribution comes 
from reaction (jSJ) or reactions ([7]) and ([6]). The overall normalization can be fit with the 
product of L p and the level of EBL. However, only the reaction ([6]) produces neutrinos. For 
a larger E max and a harder spectrum, more gamma rays come from proton pair production 
(E]), while the relative contribution required from pion producing reactions is smaller. This 
changes the flux of neutrinos relative to photons, as shown in Figs. 3 and 4. If E max is 
greater than GZK cutoff, P7cmb interactions of protons with cosmic microwave background 
radiation contribute to the neutrino flux at the highest energies. Point sources of high-energy 
neutrinos have not been obse rved so far, but IceCube detector may be able to detect such 



neutrinos (jEssey et al.ll2010bf ). 



Secondary neutrinos obey the 1/d scaling law as discussed in the introduction. This 
means that, if secondary neutrinos are detected at all, they should be detected from larger 



distances than the primary neutrinos (jEssey et al.l l2010bl ). This opens a question about 



how large the diffuse intensity from unresolved distant sources is. Obviously, it should not 
exceed direct measurements. Furthermore, considering a solid-angle bin set by the angular 
resolution of an experiment, the contribution from the diffuse background should be much 
smaller than that from a single source at a typical distance. This ensures that the apparent 
source flux is not modified and that the source can be detected as a flux excess relative to 
adjacent solid-angle bins. We find that the absolute diffuse intensity from secondaries is 
lower than that from primaries, which automatically ensures these latter requirements. 

For a simple estimate, we assume a constant density of identical, non-evolving sources 
in a Euclidean universe. For primaries, which is the familiar case, the diffuse intensity is 

dF/dQ = j-\ drLn = =7^ ( 10 ) 

J 

where n is the number density of sources (including a beaming factor as necessary), L is the 
luminosity per source, and R is the maximum distance. For secondaries, this is modified by 



including a factor (rp) in the integrand, where p is the rate of secondary neutrino production 
per unit length of the proton path: 




R 



LnR Rp 
4vr IT' 



dr Ln {rp) 



(11) 



Thus, even thought the diffuse intensity builds up more rapidly with maximum radius than 
the usual case, it is always smaller in absolute value when the probability of a proton to inter- 
act is much smaller than 1, and, therefore, Rp <C 1, the limit that we consider. Therefore, the 
questions raised above about the diffuse flux from secondaries are less important than similar 
questions for primaries, and those are addressed by the empirical definition of our model. 
In a more complete treatment, the large- R behavior would be regulated by terms describing 
evolution of sources and the expansion of the universe. For beamed sources, the beaming 
factor h should also be included; this amounts to replacing n with n/b in Eqns. f flQlfTT]) . 

One can also estimate the minimal source density n for which the diffuse flux begins to 
compete with the point sources. For an instrument with angular resolution = 1°, diffuse 
flux is given by Eqns. flUBED with dVt = 2ir (1 - cos(</>/2)) ~ 2.4 x 1(T 4 . The flux from a 
point source is L/(4nd 2 ) for primary or Lp/(4iTd) for secondary neutrinos. For a source at 
z = 0.14, or L ~ 570 Mpc, one can take R = 6 Gpc (which roughly corresponds to 2 max = 3). 
The requirement that the point source flux dominate over diffuse flux implies a condition 
on the density n of sources, or, equivalently, on the average distance I = n~ 1//3 between the 



This condition is satisfied, in both cases, for the brightest sources. 

If one can determine the flavor structure of the neutrinos from point sources, one can 
use it to further distinguish between primary and secondary neutrinos, at least in the case 
of the primary neutrinos produced in some environments with a high photon density. The 
mean free path for reaction wy — > n~p can be smaller than the decay length of the neutrons. 
Therefore, all the neutrinos would be produced from pion decays and would have the flavor 
ratios v e : : v T = 1 : 2 : 0. Neutrino oscillations alter these ratios mainly by the — > v T 
conversion, so that the detected flavors are v e : : v T — 1 : 1 : 1. However, the secondary 
neutrino flux should have a non- negligible contribution from neutron decay, n — > peu e , which 
should alter the 1:1:1 ratio. Of course, the observed signal could be a combination of both 
primary and secondary neutrinos. The flavor combination of such a combination should still 
differ from 1:1:1, but the spectra include an additional component. Thus, we expect the 
signals plotted below to be the lowest neutrino signal to be expected. 



sources: 




14 (10 3 /6) 1 Mpc, for primary neutrinos 

1 /3 

8 (10 3 /6) Mpc, for secondary neutrinos. 



(12) 
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Fig. 3. — Gamma-ray spectra (lower-energy curves) and neutrino spectra (higher-energy 
curves) for different values of maximal proton energy, -E max = 10 8 GeV, 10 10 GeV, 10 11 GeV 
(top to bottom), for sources at red shifts z = 0.14 (left), such as blazar 1ES 0229+200, and 
z = 0.44 (right), s uch as blazar 3C66A. T he cosmic ray luminosity L p was adjusted to fit th e 



data from HES S (lAharonian et al 



and VERITAS flAcciari et al. 



2008 



2006h. MAGIC JAlbert et~al1l2008l : 



Aliu et al. 2009 



Acciari et al 



2009J) 



Aleksic et all 120101 ) 



Individual curves 

are labeled by the value of the spectral index a. Here we assume vanishing IGMFs; a ~fG 
or higher magnetic field would cause some reduction of flux below 1 TeV. 
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Fig. 



4. — Gamma-ray spectra (lower-energy curves) and neutrino spectra (higher-energy 
curves) for a = 2.0,2.3, and 2.6 (top to bottom), for sources at red shifts z = 0.14 (left) 
and z = 0.44 (right). The cosmic ray luminosity L p was adjusted to fit th e data from 
HESS jAharonian et all 120061). MAGIC ([Albert et al.ll2008l: kleksic et al.l fioioh. and VERI- 



TAS (Acciari et al. 



2008 



Aliu et al.ll2009l ; lAcciari et al.ll2009l ). Individual curves are labeled 
by the value of the maximal proton energy -E max - Here we assume vanishing IGMFs; a ~fG 
or higher magnetic field would cause some reduction of flux below 1 TeV. 
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Implications for extragalactic background light and cosmic ray acceleration 

models 



Our predictions for the gamma-ray spectra fit the data extremely well for all the models 
of EBL. While primary gamma rays are lost to interactions with EBL, secondary gamma 
rays are produced in these interactions. The overall flux depends on the EBL photon density 
multiplied by the unknown and poorly constrained L p , the energy output of AGN in protons. 
There are small differences in spectral shapes, and one can hope to gain some discriminating 
power with more data. At present, the data shows some preference for high EBL, although 
it is not statistically significant (see Table 1). (This is in contrast with the limits set on 
EBL under the assumption that all the observed gamma rays are primary.) The inclusion of 
secondary gamma rays brings in one additional "free" parameter L p , but it affects only the 
overall normalization of the spectrum. The shape of the spectrum fits the data quite well, 
much better than the fit one could obtain with the primary sources alone. Moreover, fitting 



be oossible (Stecker et al. 


200' 


7 ), but which are by no means natural or generic, based on a 


number of theoretical models ( . 


^eacock 


1981 




£irk & Schneider 1987; Heavens & Drurv 1988; 


Bednarz & Ostrowski 1998 


; Malkov & Drurv 


2001 


). The data thus favor the interpretation 



in terms of secondary gamma rays, which relaxes the constraints one can put on the models 
of EBL. Cosmic ray acceleration models can be improved if the neutrino observations provide 
insights about the maximal energy and spectral slope. 

In general, secondary gamma-ray flux depends on the present level of EBL more than 
the EBL evolution, while the neutrino fluxes probe both present and past levels of EBL. The 
detection of point sources by IceCube, combined with improved gamma-ray data, can help 
improve the bounds on both the present level of the EBL and on its evolution. 

In principle, one can use gamma-ray data to set upper and lower limits on cosmic ray 
production in AGN by comparing the primary gamma-ray component with the component 
generated by cosmic-ray interactions along the line of sight. There is no doubt that the 
former dominates the signals from nearby blazars, while the latter can take over at large 
distances, where the former is filtered out by photon-photon interactions. However, since 
the predicted gamma-ray spectra show little sensitivity to the values of the lower and upper 
energy cutoffs and the spectral slope for cosmic rays, the ratio of the power in cosmic rays to 
the power in gamma rays can vary dramatically, depending on these uncertain parameters. 
We have not been able to set meaningful limits based on the present data. Uncertain as they 
are, the li mits on the power in c osmic rays derived from comparison of models with cosmic 
ray data (IBerezinsky et all 120061 ) appear to be more constraining than any limits one could 
derive from the present gamma-ray data. 
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Table 1. Model parameters for the spectra shown in Fig. 1. (Here we assumed 
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Fig. 5. — Dependence of the isotropic equivalent of source power L p> iso in cosmic rays on the 
proton spectral index a and the maximal proton energy E max for sources at two redshifts, 
1ES0229+200 at z = 0.14 (left) and 3C 66A at z = 0.44 (right). Here IGMF effects are 
neglected (which results in a small difference with the values given in Table 1). 
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Conclusion 



The surprisingly low attenuation of high-energy gamma rays can be explained by sec- 
ondary gamma rays produced in interactions of cosmic-ray protons with background photons 
in the intergalactic space. We have obtained excellent fi ts to observed spectra of several dis- 



tant b lazars, hence extending our prior published results (jEssey &: Kusenkoll2010l ; lEssey et al. 



2010b|). All these spectra have a characteristic shape that derives from the shape of the EBL, 
and they are less sensitive to the level of EBL than the spectra of primary gamma rays. At 
low energies the spectra a re harder than predicted by theoretical models ( jMalkov fc Drury 
200ll ; IStecker et al.l 120071 ) , which explains why Fermi has not detected these high-energy 
sources. A future detection at low energies could help differentiate between the primary and 
secondary gamma rays. Secondary gamma rays are expected to show no temporal variability, 
which so far is consistent with the data at energies above TeV. The temporal information is a 
strong discriminant between primary and secondary gamma rays; most primary gamma-ray 
models predict temporal variability on a time scale smaller than the observation time. 



We have also presented our predictions for secondary neutrino signals from blazars. The 
power for the neutrino signal peaks at about 1 — 10 eV cm~ 2 s _1 , depending on the choice 
of source and model parameters. For the declinations of the sources we considered, IceCube 
sensitivity is, roughly, 1 eV cm~ 2 s -1 for 2 2 strings after 0.75 years and 2 eV cm _2 s _1 for 



80 strings after 1 year (lAbbasi et al.l 120091 ) . This makes the prospects for detection seem 



plausible. However, the IceCube sensitivity is calculated assuming a standard E~ 2 spectrum 
in the energy range from 3 TeV to 3 PeV. A detailed analysis of IceCube sensitivity for the 
specific spectrum shape predicted for secondary neutrinos and the higher energy range is 
beyond the scope of this paper. 

One can make additional predictions for neutrinos, besides the spectral shape. First, 
as for gamma rays, there should be no temporal variability observed for neutrinos. Second, 
the luminosity of sources should vary with distance as 1/d, as opposed to the usual 1/d 2 
scaling law. This makes more distan t sources observable, as compared to predictions for 
primary neutrinos (IStecker et al.lll99ll ). and it allows one to confirm the secondary neutrino 
observations by studying a population of distant sources. Finally, the flavor structure of the 
observed signal should differ from primary neutrino models due to a significant contribution 
from neutron decays. 

Secondary gamma rays and neutrinos present a new powerful method to probe the 
radiation and magnetic field contents of intergalactic space, as well as AGN properties. The 
signals differ from primary signals in spectral shape, temporal variability, and scaling with 
distance. Future ACT and neutrino experiments should consider the effect of specific scaling 
laws for secondary gamma rays and neutrinos, which bring more sources into the field of 



view of a given instrument. 
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